Metabolic and physiological effects of water stress on Moshgak (Ducrosia anethifolia Boiss) populations using GC–MS and multivariate analyses

Ducrosia anethifolia Boiss is a medicinal and aromatic plant that is used as a food additive and drug in the areas of natural distribution. In this study, eight populations from different geographical regions of Iran were evaluated for their essential oil composition, total phenolics and flavonoids as well as for free radical scavenging activity during two consecutive years under water deficit. Analysis of variance was performed using a randomized block design at three levels of irrigation and 2 years, using the GLM procedure of SAS software and cluster analysis was used according to Ward's method using squared Euclidean dissimilarity. The highest essential oil and total phenolics/flavonoids content were obtained in the second year. GC–MS analysis was used to analyze the essential oil components. In normal conditions, cis-chrysanthenyl acetate (39.1–66.4%) and α-pinene (1.02–4.7%) were recognized as major compounds. The essential oil components were elevated under water stress. The highest increase in cis-chrysanthenyl acetate (21.32%) in response to water stress was observed in Abarkuh1. Elevation in total phenolic, flavonoid and antioxidant activity was also observed in plants subjected to water deficit. The highest content of phenolic acids was measured in Semirom (136.74 mg TAE/g DW), while the highest flavonoid content was in Kerman (6.05 mg QUE/g DW) in severe water stress condition. Finally, a high diversity in the studied populations can be used to select advantageous populations for other pharmaceutical and food purposes.

www.nature.com/scientificreports/ Today, the uses of many bioactive components such as phenolics and flavonoids is of the great value for human health 6 . Recently, many studies have reported the composition of essential oils, phenolics and flavonoids as well as their antioxidant activities in different medical herbs. The number of compounds and their bioactive capacities might be affected by genetic as well as environmental factors 8,9 . Drought is one of the main environmental factors that affects the concentration of metabolites 10 . Drought or osmotic stress leads to excessive generation of reactive oxygen species (ROS) that can be harmful to cellular components at high concentrations 11 . Plants demonstrate a wide range of stress tolerance or stress avoidance mechanisms to survive in a stress condition 10 . Plants under drought stress conditions mostly produce higher content of secondary metabolites compared with non-stressed condition 12 . There are many studies that highlight the effect of water stress on secondary metabolites of Apiaceae plants, including Cuminum cyminum 13 and Foeniculum vulgare 14 . Total phenolic (TPC) and total flavonoids (TFC) and antioxidant activity can also be elevated in response to water stress in different medicinal plants 1,8,14,15 . Also, Mottaghipisheh et al. 16 showed that the antioxidant activity of D. anethifolia in ethanolic and ethyl acetate extracts is less than Butylated hydroxytoluene (BHT) as a synthetic antioxidant that was used as positive control. Evaluation of populations in various years is one of the substantial segments of the breeding programs. Moreover, cultivation year can highly affect the secondary metabolites of medicinal plants such as Foeniculum vulgare Mill 9,17,18 and Stachys lavandulifolia 19 .
There are limited studies on D. anethifolia populations. Sabbaghi and Rahimmalek 20 assessed the morphological and molecular variation of Iranian Moshgak populations. Arbabi et al. 2 also evaluated morphological and oil components of leaves in this species in one geographical region in Southeastern parts of Iran. However, no studies were performed on the effect of water stress on the composition of essential oils, phenolics, flavonoids, and antioxidant capacity of different populations of D. anethifolia in two consecutive years. As gas chromatography results provide valuable data with many information that is difficult to process based on univariate methods, multivariate analyses provide new tools for exploration and interpretation of chromatographic data [21][22][23] .
The purposes of this research were: (1) to evaluate different populations of D.anethifolia under experimental field conditions to select populations with high phytochemical contents; (2) to assess chemical composition and polyphenolic content in water stress condition and (3) to use multivariate analyzes for choosing and introduction of elite populations for further cultivation and processing.

Results
Effects of water stress and population on essential oil content. The ANOVA results showed significant differences between the Moshgak population studied for the essential oil content. Furthermore, the effect of year, irrigation, year × population interactions and irrigation × population interactions were significant for EOs ( Table 1). The color of the extracted essential oil content was light yellow. However, the trends of essential oil over two consecutive years ( Table 2) showed that the higher amount of the EOs occurred in the second year. The EOs content also increased under drought stress (Table 1).
In the first agronomic year (2019), the highest and lowest increase in essential oil content was observed for populations of Abadeh (83%) and Ferydunshahr (52%), respectively (Table 2).
In the second agronomic year (2020), the water deficit caused significant increases in essential oil in all Moshgak populations. The highest EOs enhancement was observed in the Qazvin and Abarkuh2 populations (83%), while the lowest in Abarkuh1 (56%) ( Table 2).
Effects of drought stress and population on essential oil composition. The essential oil constituents are listed in Table 3. Chromatogram of GC-MS peaks corresponding to D. anethifolia and the mass spectra of major compounds are illustrated in Fig. S1. The results indicated a high chemical polymorphism among Moshgak populations, and the 26 EOs components were affected both by irrigation regime and population of origin. The main constituents detected in the EOs were cis-chrysanthenyl acetate and α-pinene, ( Table 3). The results also showed that cis-chrysanthenyl acetate was the most abundant component. In control conditions, the highest and lowest amount of this compound was observed in Semirom (66.4%) and Abadeh (39.1%), respectively. Furthermore, the concentrations of the essential oil composition in all Moshgak populations elevated Table 1. Analysis of variance (mean squares) for phytochemical traits of 24 Moshgak populations (P) evaluated at three levels of irrigation (I) and 2 years (Y) of 2019 and 2020. DF: degrees of freedom; ns, nonsignificant; *, significant at P ≤ 0.05; **, significant at P ≤ 0.01.  Table 3). The highest and lowest increases in cis-chrysanthenyl acetate were in Abarkuh1 (21.32%) and Hormozgan (6.42%), respectively (Table 3). Under water deficit, α-pinene the highest increase was in Kerman populations (6.1%) and the lowest in Qazvin (0.61%) ( Table 3). The maximum increase in trans-verbenol was observed for the Abadeh population. Sabinene was not present in the Abadeh, Ferydun-shahr1, and Abarkuh1 populations. However, in fruits from drought stressed plants the emergence of sabinene noticed ( Table 3).

Effects of drought stress and population on phenolic contents.
The results of the analysis of variance revealed a significant variation between the Moshgak population studied for the total phenolic content (TPC) and flavonoid content (TFC). In the present study, TPC and TFC were affected by year and irrigation ( Table 1). The interaction effects were significant between the population × irrigation and the year × irrigation × population for these characteristics. During the first season of experiment, under control conditions, "Ferydunshahr" exhibited the highest TPC (95.87 mg of TAE/g of DW), while 'Kerman' showed the highest TFC (4.13 mg of QE/g of DW). Elevated total phenolic and flavonoid content in all populations under drought stress (Table 2).
In the second year, under control conditions, the highest TPC was observed in Ferydunshahr (99.34 mg of TAE/g of DW), whereas the maximum TFC was recorded in Kerman (4.68 mg of QE/g of DW). Drought stressed samples from all populations had significantly increased total phenolic and flavonoid content. The highest increase of TPC and TFC were in populations of Abarkuh (52.87%) and Kerman (1.37%), respectively ( Table 2).
Effects of drought stress and population on DPPH radical scavenging assay. The results of ANOVA indicated a significant difference among the Moshgak population studied for the DPPH radical scavenging assay (IC 50 ). Furthermore, the effect of year, irrigation, population year × population irrigation × population interactions was significant for IC 50 (Table 1). In the first year, under control conditions, the highest antioxidant activity (the lowest IC 50 value) was obtained in Hormozgan (166.17 µg/ml), while Abarkuh1 possessed the lowest antioxidant activity (the lowest IC 50 value) (726.81 µg/ml). Under drought stress antioxidant activity increased (decreased IC 50 ) in all Moshgak populations ( Table 2). Table 2. Mean comparisons of irrigation regimes, years and populations (interaction) on phytochemical traits of 8 Moshgak population. In each trait, means having difference less than the LSD provided are not significantly different (p < 0.05). www.nature.com/scientificreports/ In the second year, under control conditions, the lowest and highest IC 50 were observed in Hormozgan (216.63 µg/ml) and Semirom (730.13 µg/ml), respectively (Table 3). Additionally, the activity increased in response to drought stress ( Table 2).

Irrigation regimes Population
Correlations between the main essential oil components and TPC and TFC. The correlation coefficients among the main components of essential oil and total content of phenolics and flavonoids are presented in Table 6. In control conditions, the highest coefficeints were between α-pinene and β-pinene (r = 90**), as well as myrcane and limonene (r = 84**). Cis-chrysanthenyl acetate had a weak positive correlation with β-pinene (r = 0.66*), but was slightly negatively correlated with limonene (r = − 0.46*) ( Table 6). Furthermore, the relationship between p-cymene and total phenolics (r = 66**), and IC 50 (r = 0.48 *) was positive, but the correlation was negative with essential oil content (r = − 0.54*). However, there was no significant correlation between cischrysanthenyl acetate and flavonoids or phenolics content (Table 4).
Hierarchical cluster analysis and principal component analysis. To   www.nature.com/scientificreports/ the average of 2 years (Fig. 1). The HCA results can be applied to choose the superior parents based on the suitable genetic distance and cis-chrysanthenyl acetate for targeted crossing in future Moshgak breeding programs. According to the HCA, the Moshgak populations studied were grouped into three groups. Cluster I comprised 12 populations of normal and moderate water stress conditions. The second group included Abarkuh populations from moderate and severe water stress conditions with the highest cis-chrysanthenyl acetate, p-cymene, geranyl acetate, total phenol and IC 50 . Eight populations from moderate and severe water stress conditions were placed in the third group. They were distinguished by high α-pinene, sabinene, β-pinene, myrcane, limonene, γ-terpinene, trans-verbenol, decanal, and total flavonoids (Fig. 1). The principal component analysis (PCA) was performed based on the main phytochemical traits averaged for 2 years. The principal component analysis confirmed the results of HCA (Fig. 2). The PCA result indicated that the first two components explained 39% and 19% of the total variation, respectively (Fig. 2). The first PC (PC1) had a positive association with α-pinene, sabinene, β-pinene, myrcane, limonene, γ-terpinene, transverbenol, decanal and flavonoids. The second PC (PC2) had a high correlation with the essential oil content, cis-chrysanthenyl acetate, p-cymene, geranyl acetate and total phenolic and IC 50 . The Abarkuh1 and Abarkuh2 populations in moderate and severe water stress conditions were classified in the same group. Abadeh and Kerman populations under moderate and severe water stress conditions were classified near each other. Additionally, Ferydunshahr, Semirom, Hormozgan,, and Qazvin populations under severe water stress conditions were classified in the same group. The others are classified into a separate group. Based on the biplot analysis, selection according to the high PC1 and PC2 values leads to populations with high cis-chrysanthenyl acetate, geranyl acetate, and total phenolics (Fig. 2).

Discussion
Medicinal herbs have been used for thousands of years by many people all over the world. Due to the overharvesting of natural habitats, the expansion of procedures for preservation, assessment and application of medicinal herbs is essential for the protection and utilization of genetic resources 18 . Narrow information on genetic diversity and the effect of drought stress on Moshgak is available. In this study, we analyze various Moshgak populations collected from different geographical regions of Iran by phytochemical traits and essential oil components under drought stress conditions. The results of this study presented significant differences among populations and irrigation regimes for all characteristics studied, representing a wide diversity among populations in response to drought stress in this germplasm. Previously, a wide variation in important morphological traits and oil content of local Moshgak populations was reported 2,7 . The results of this research showed that populations are significantly different for EOs and reacted differently in response to drought stress. The increase in EOs production in response to a water deficit may be due to carbon reabsorption, as plant growth is reduced, for the biosynthesis of stress-related metabolites 24 . Previous studies reported different ranges of essential oil content in leaves of Moshgak populations 16,25 , while in the present study, the oil content ranged from 1.31 to 3.68% in fruits. EOs content can be influenced by environmental factors 4,19 . Alinian and Razmjoo 26 observed that drought stress can lead to an increase in Cumin EO (Cuminum cyminum L.). In addition, Bettaieb et al. 13 reported that fruits essential oil decreased in cumin in response to severe drought stress, while it increased in response to moderate conditions of water deficit. These findings indicated that the EOs depend on genotypes, and drought level. The results of the GC-MS analysis indicated a wide chemical polymorphism among the Moshgak populations. In the present study, the oil components fruits of this plant were determined for the first time. Arbabi et al. 2 reported Table 4. Pearson correlation coefficients between major essential oil compounds of fruit essential oil, total phenol, total flavnoid and DPPH of eight Moshgak population in control (lower triangle) and moderate water stress (upper triangle) conditions. *Significant at P ≤ 0.05; **significant at P ≤ 0.01. www.nature.com/scientificreports/ Table 5. Pearson correlation coefficients between major essential oil compounds of fruit essential oil, total phenol, total flavonoid and DPPH of eight Moshgak population severe water stress conditions. *Significant at P ≤ 0.05; **significant at P ≤ 0.01. www.nature.com/scientificreports/ that cis-chrysanthyl acetate, n-decanal, decanol, and α-pinene were the main compounds of leaves in populations of D. anethifolia. Furthermore, Karami and Bahloli 27 reported n-decanal, dodecanal, cis-chrysanthenyl acetate, and α-pinene as the main compounds in the essential oil of D. anethifolia leaves. In the present study, cis-chrysanthyl acetate and α-pinene were the main compounds. It seems that the production of cis-chrysanthenyl acetate starts only after flower appearance. As the production of cis-chrysanthenyl acetate in the leaves was very low, the biosynthesis of this compound was carried out mainly in the full grown inflorescence 28 . Prostaglandin synthetase is inhibited by chrysanthenyl acetate and could have analgesic properties 29 . Furthermore, chrysanthenyl acetate might contribute to antimigraine activity, due to its inhibition of prostaglandin 30 . Cis-chrysanthyl acetate is a monoterpene hydrocarbon 2 . Previous studies revealed that monoterpenes can be highly influenced by environmental stress 31 . Plants produce high concentrations of monoterpenes under environmental stress conditions because of a low allocation of carbon to growth, suggesting a competition between defense and growth 32 .
In the present study, some of essential oil compounds increased in response to drought stress. In line with our findings, Bettaieb et al. 13 reported that the application of irrigation regimes can cause significant changes in the cumin EOs quality. For example, the percentage of γ-terpinene increased significantly in response to drought stress. The variation in EO compounds at various levels of water deficits is due to changes in the activity of linked biosynthesis enzymes in response to water deficit 13 . Monoterpene formation is catalyzed by terpene synthases, which are activated by developmental and stress-linked factors 33 .
In this study, the total phenolic and total flavonoid content increased in response to water shortages, though they pursued various trends, suggesting that the biochemical responses of Moshgak to drought stress could be population dependent. Similarly, the phenolic contents of cumin under water stress increased greatly 13 . In contrast, Neffati et al. 34 reported that the polyphenol content decreased under salinity stress in coriander (C. sativum L). The different trends in phenolic fractions under drought stress can be imputed to the synthesis of various components of TPC and TFC as efforts by various populations to adapt to different soil moisture levels. Many researchers reported that the accumulated components of TPC and TFC could be increased under drought  www.nature.com/scientificreports/ stress as a means of improving crop quality in medicinal plants 7 . Harrison and Were 35 reported that improved TPC and TFC in response to drought stress could be the result of the degradation of bigger phenolic components into smaller ones. Furthermore, phenolics have strong antioxidant activity, and therefore participate in advocacy versus reactive oxygen species (ROS) that are inevitably produced when photosynthetic metabolism is devastated by environmental constraints 13 . The results of the present study indicated that the populations are conspicuously different for IC 50 and reacted differently in response to drought stress. Furthermore, variation was observed in the antioxidant activity of drought-stressed plants according to those of their TPC. In corroborating this finding, Bettaieb et al. 36 indicated that the variation in the antioxidant activity of cumin was reflective of the variation in its TPC. Other researchers reported that cumin antioxidant activities exhibited that treated seeds with moderate and severe water deficit showed the maximum antioxidant activity of cumin 13 . Exposure of plants to adverse environmental conditions like salinity, extreme drought, or temperature stresses can enhance the production of ROS 13 . Low molecular weights antioxidants can effectively scavenge harmful radicals 11 . In the present study, positive relationships were also observed between IC 50 and TPC under drought conditions. Compared to our finding, previous studies indicated a significant correlation between antioxidant activity and TPC in C. cyminum 13 and C. sativum 34 .
The results of the present study led to identified populations with high cis-chrysanthenyl acetate, essential oil from seeds, and total phenolic content using PCA (high PC1 and PC2). The Abarkuh populations under stress were classified into a separate group due to the high p-cymene content. Hosseini et al. 18 reported that selection based on a biplot led to the identification of genotypes with high essential oil yield and main components in fennel.
Various relationships were obtained among the characteristics. Positive correlation of some compounds showed that selection according to one of the compounds can cause improvement of another one 18 . The results of the present study revealed that α-pinene had a negative correlation with cis-chrysanthenyl acetate. In addition, a negative correlation was observed between cis-chrysanthenyl acetate with trans-verbenol and limonene. Negative associations between some isomer compounds were also reported in previous research 4, 6 .
The water deficit induces partial stomatal closure and massively reduces the CO 2 influx into the leaves, which directly leads to a passive increase of all processes consuming NADPH + H + , including the biosynthesis of highly reduced secondary plant products. Alternatively, enzymes responsible for the biosynthesis of natural products could be actively up-regulated. The corresponding increase in NADPH + H + consumption significantly leads to the dissipation of the stress-related energy and thus prevents the generation of oxygen radicals. These coherences reveal the stress mechanism that leads to increased biosynthesis of natural products such as isoprenoids, alkaloids, or phenolics 37 .

Conclusion
The Moshgak germplasm used in the current study indicated a wide range of variations in phytochemical traits and essential oil compounds that can be useful for further breeding programs. The total phenolic and total flavonoid content, antioxidant activity and EOs increased in response to water shortages. Also, the essential oil compounds increased in response to drought stress. Furthermore, cis-chrysanthenyl acetate was the most abundant essential oil component in the fruits of this plant. The highest increase in cis-chrysanthenyl acetate (21.32%) in response to water stress was observed in Abarkuh1.Significant differences in response to drought stress revealed that each population can be exploited for certain purposes.
Superior populations introduced in the current research are suitable candidates for the further development of new cultivars for medicinal and pharmaceutical programs.

Materials and methods
Ethics statement. This study was carried out in accordance with the Ethics Committee of Isfahan University of Technology, Iran. Written, informed consent was obtained from all volunteers. All experimental protocols were approved by the Ethical Committee in Isfahan University of Technology, Iran.  Table 6 and Fig. 3. The experiment was carried out on the IUT Research Farm, Esfahan (Latitude of 32° 38ʹ North, Longitude of 51° 39ʹ East). Seeds of all populations were planted in a Randomized Complete Block Design (RCBD with three replicates and three irrigation regimes (three treatments). Each plot included two rows of 2 m in length and 70 cm apart. The soil was clay loam texture, pH 7.5 and organic matter content of 1%. Fertilizers were used at 100 kg N/ha and 100 kg P/ha prior to sowing. Irrigation regimes were specified according to the procedure of Allen et al. 37 . The experimental research and field studies were performed based on the standard methods as described by the previous literatures 18,32,36 . Water stress treatment. Before irrigation, soil samples were obtained at soil depth of 0-15, 15-30, and 30-45 cm from three environments to estimate the soil water content and to compute the irrigation water content based on a 45-cm rooting depth. The irrigation depth was estimated according to the following formula:  38 . Seeds were sown under three experimental conditions, including control condition (C), moderate (MS), and severe water stress (SS). The plants were irrigated at the sowing time, after which irrigation was carried out every week until the early flowering stage (20% flowering). Moreover, irrigation was continued until 40, 60 and 80%, respectively, of the field capacity (FC). To obtain the volume of each water application, a Parshall flume formula (Id = I × p; where p is the fraction of I that may be depleted from the root zone) was applied. Finally, Ig = (Id/Ea) × 100, with Ea representing irrigation efficiency (%), was supposed to be 65% on average.
Essential oil extraction. Ripe fruits from D. anethifolia plants grown in each agronomic year were collected and dried in the shade at 25 °C for 4 days. The fruits of ten samples from each population were mixed and then applied to measure essential oils. The oil was obtained by hydro-distillation using a Clevenger-type apparatus for 4 h. For this purpose, 50 g of each samples was applied. Distilled essential oil was procured using diethyl ether as collecting solvent (v/v), collected in the glass container, and kept at 4 °C until analysis. Finally, the essential oil content was presented as percentage based on the weight of the air-dried seed specimens.
Gas chromatography-mass spectrometry (GC-MS) analysis. Essential oil compositions were analyzed using an Agilent 7890 gas chromatograph (Agilent Technologies, Palo Alto, CA, USA) with a HP-5 MS capillary column (30 m × 0.25 mm, 0.25 μm film thickness). The oven temperature was programmed at 70 °C for 2 min followed by an increase to 200 °C (Hold = 1 min). Injector and detector temperatures were set at 210 and 220 °C, respectively. The analyzes were done via helium as the carrier gas (inlet pressure 17.7 psi; split 1:10). The gas chromatograph was coupled to the Agilent 5975C mass spectrometer as a selective detector with an ionization voltage of 70 eV; ion source and mass quadrupole temperatures of 230, and 150 °C, and mass range of 35-500 m/z. The GC MSD Chem Station program was applied as operating software. Retention indices were obtained by using retention times of n-alkanes (C8-C21) injected before essential oil analysis in the same conditions, using the Kovats index formula 39 . Essential oil constituents were identified by comparison of retention indices (RIs) with those reported in the literature 39 , as well as comparing their mass spectra with those recorded in the NIST 08 (National Institute of Standards and Technology), Wiley275. L (ChemStation data system). The percentage of composition was obtained on the peak areas of the total oil compositions.
Methanolic extract and total phenolic content. The total phenolic content was obtained with the Folin-Ciocalteu reagent as explained in Hodaei et al. 8 with some changes. Briefly, to provide a sample extract, 25 mL of 80% methanol was added and shaken at 25 °C for 24 h. The extract was filtered, and 2.5 mL of the Folin-Ciocalteu reagent (1:10 diluted with distilled water) and 2 mL of 7.5% sodium carbonate was added to 0.5 mL of the methanolic extract and shaken well. The mixture was heated at 45 °C for 15 min. Finally, the absorbance of the extract was read at 765 nm via spectrophotometer. The water and reagents was applied as a blank sample. Tannic acid was applied as the reference standard and the total phenolic content was reported as mg of tannic acid equivalents (TAE) equivalent per gram of sample on a dry basis.  www.nature.com/scientificreports/ Total flavonoid content of the extracts. The aluminum chloride colorimetric procedure was used for the determination of total flavonoids by Sarfaraz et al. 40 with minor changes. The methanolic extract (125 µL) was mixed with 75 µL of a 5% NaNO 2 solution and after six minutes of keping in the dark, 150 µL of AlCl 3 (10%) was added and incubated for 5 min. Then, 750 μL of NaOH (1 M) and 2500 μL of distilled water were added. After 15 min of incubation, the extract turned pink and the absorbance was read at 510 nm. Finally, TFC was reported in mg of quercetin equivalents (QE) per gram of fruits dry weight.
DPPH radical scavenging assay. The radical scavenging activity of D. anethifolia was determined using the procedure of by Sarfaraz et al. 40 with minor changes. Briefly, 5 mL of 0.1 mM methanol DPPH solution with 0.1 mL of the extract was mixed at selected concentrations (50, 100, and 300 ppm). The mixture was shaken and kept for 30 min at 25 °C. The absorbance of the samples was read using a UV-visible spectrophotometer at 517 nm with methanol (80%) as a blank and butylated hdroxytoluene (BHT) applied as standard antioxidant at 1-100 μg/ml for comparison. The absorbance was measured at 517 nm versus the blank. The percentage of inhibition of the samples was estimated on the basis of the following formula: where, A A and A B are the absorbance values of the DPPH radical in percentage in the fruits extract sample and the control, respectively. Subsequently, the inhibition percentage was plotted against the extract concentration and 50% of the inhibitory concentration (IC 50 ) of the DPPH values was calculated (Fig. S1).

Statistical analysis.
The normality of the numerical data was computed using the proc normality test of SAS software release 9.2 (SAS Institute, Cary, NC, USA) 41 . Analysis of variance (ANOVA) was performed using a randomized block design at three levels of irrigation and 2 years, using the GLM procedure of SAS software release 9.2 (SAS Institute, Cary, NC, USA) 37 . Mean comparisons were made using the least significant difference (LSD) test (p < 0.05) using SAS software release 9.2 (SAS Institute, Cary, NC, USA) 37 .
Cluster analysis was used according to Ward's method using squared Euclidean dissimilarity through the Stat Graphics software ver. 17.2 (http:// www. statg raphi cs. com). Principal component analysis (PCA) was performed to recognize the association among D. anethifolia populations and all measured characteristics using Stat Graphics software ver. 17.2 (http:// www. statg raphi cs. com).

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author based on reasonable request.